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Abstract 


Impedance matching networks have wide range of application in the field 
of communication. Recent survey of literatures show different design 
procedures to realize its design and enhance its performance. This paper 
evaluated four nature-inspired optimization algorithm and introduces a 
new optimization algorithm, Modified Firefly Algorithm (MFA), which 
improves performances of the Standard Firefly Algorithm (SFA) in many 
ways. MFA was used to design both lossless and lossy impedance 
matching network and yielded results that were comparable to other 
optimization algorithms in literatures. 


Index Terms: Global Optimization Method, Impedance Matching Network, 
Firefly Algorithm. 
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Antenna Impedance Matching 


m [mpedance matching is the connection of additional 
impedance elements to existing ones in order to achieve 
maximum power transfer, or to reduce reflection in the 
transmission line. 


Wire Dipole Antenna cai 


impedance 


Matching Transmission Line 
System 


Radio Receiver 





Typical Impedance Matching System 
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Performance of Impedance Matching 
Netwotk (IMN) 


m The performance of the IMN can be evaluated 
by a parameter called Transducer power gain 
(TPG). 

m The other goal is to minimize the reflection 
coefficient (/'). We could use the VSWR for 
this objective since it is directly related to the 
reflection. 
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Impedance Matching Technigues 


Impedance Matching 
for Antenna 





Distributed 
Structural Modification 


Lumped element 


CAD Method 






Classical Nature-Inspired 


Impedance matching Technique diagram 
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Outline 


m |ntroduction to impedance matching 

m Nature-inspired optimization algorithms 
m Modified firefly algorithm (MFA) 

m MFA on impedance matching problem 


Lossless impedance matching 
Lossy impedance matching 
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Optimization Algorithms 


m Classical Methods 
Newton’s Method 
Gradient Descent Method 


m Nature-Inspired Methods 
Genetic Algorithm (GA) 
Particle Swarm Optimization(PSO) 
Ant Colony Optimization (ACO) 
Firefly Algorithm (FA) 
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Genetic Algorithm (GA) 


m Concept of GA 


(1,2) (1,3) (1,4) (1,5) 











, (3,4) 
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Swatm Intelligence Algorithm 
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Swarm Intelligence Algorithm 


= Swarm intelligence algorithms are derived from the 
social behavior of animals that consist of a group of 
non-intelligent simple agents with no self-control to 
their behavior. The local interaction between the 
agents and their environment lead to the swarm 
behavior that is new to the individual agents. 


m Examples of swarm intelligence algorithms are PSO, 
ACO, and FA. 
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Particle Swarm Optimization (PSO) 


m Concept of PSO 





Iteration # O Iteration # N 


Vi = Vi + cı * rand * ( Dpest - Xj) + C2 * rand * (gpest — Xi) 


Xj; = Xj + Vi 
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Ant Colony Optimization (ACO) 
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Firefly Algorithm (FA) 
m Attraction of a firefly to another is 
proportional to the Light intensity. 
Light intensity 


Is 
I(r) = 2 
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Lossless Impedance Matching 


m The Pi3 network is shown below 


Antenna 
Z 





Lossless Pi3 network 


m This network was used for an impedance matching network. The 
load impedance can be obtained from measurement of a 80-meter 
inverted-V dipole antenna between frequencies 3.5-3.85 MHz. In 
simulation the antenna load data is from a record and a source 
resistance is 50 ohms. 
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Lossless Impedance Matching 


m 7ransducer Power Gain 
TPG = 1 — |r|? 


m J’ isthe reflection coefficient, such that 
l = (Zin — Zs )/( Zin + Zs) 
1+r 
1-r 
where Z;, is the input impedance seen by the 
generator at the equalizer input terminal. 


VSWR = 
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Optimizers 


m Global Optimization Toolbox from the MATLAB 
Simulated Annealing (SA) 
Pattern Search (PS) 
Genetic Algorithm (GA) 


m MATLAB codes of Nature-lnspired Optimization Methods 
Genetic Algorithm (GA) 
"| Particle Swarm Optimization (PSO) 
Ant Colony Optimization (ACO) 
Firefly Algorithm (FA) 
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Results — Global Optimization Ioolbox 


VSWR with matching network optimized 












VSWR with matching network optimized 45 
—— 100 loops 
—— 100 | 
= =p 
oops 4 
35 
3 
ć 
2.5 
2 
15 
1.4 L 1 L 1 4 
JA as a8 ad a8 a2 Š Mu 35 38 37 38 38 4 
Freguency/Hz x 107 Frequency/Hz x 107 
Pattern search Simulated Annealing 
VSWR with matching network optimized VSWR with matching network optimized 
35 
—— 100 loops —— 100 loops 
—— 500 loops —— 500 loops 







—— 1000 loops 
——— 5000 loops 





ba 35 3.6 37 3.8 33 4 
Frequency/Hz x 10° 


GA with population = 20 





—— 1000 loops 
——— 5000 loops 


35 3.6 37 3.8 AS 4 
Frequency/Hz x 107 


GA with population = 200 
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esults — Nature-Inspired Methods 


VSWR with matching network optimized VSWR with matching network optimized 








1 
ha 3.5 3.6 37 3.8 3.9 4 La 3.5 3.6 3.7 3.8 3.9 4 
Frequency/Hz x 10? Frequency/Hz x 10° 

VSWR with matching network optimized VSWR with matching network optimized 





14 3.5 3.6 3.7 3.8 3.9 
Frequency/Hz x10 Frequency/Hz 


ACO FA 


La 3.5 3.6 3.7 3.8 3.9 4 
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Results - Comparison 


BT GLOBAL OPTIMIZATION TOOLBOX Nature-Inspired Optimization Methods 
Algorithm PS SA GA 


N/A N/A 20 200 70 150 70 150 70 150 70 150 
100 100 100 100 N/A N/A N/A N/A N/A N/A N/A N/A 


1.4758 1.8664 2.0386 1.9199 1.3386 1.3386 1.3386 1.3386 1.3386 1.3386 1.3386 1.3386 
VSWR 


1% error 0 0 0 0 981 995 512 807 521 833 905 907 







VSWR<1.352 
2 %error 0 0 0 0 997 1000 618 860 631 884 985 983 


VSWR<1.365 
3 % error 0 0 0 0 1000 1000 984 990 692 901 998 998 
VSWR<1.379 


11.00 154.00 29.00 186.00 1.16 1.11 1.11 1.17 1.16 1.14 1.14 1.20 
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Outline 


m |ntroduction to impedance matching 

m Nature-inspired optimization algorithms 
m Modified firefly algorithm (MFA) 

m MFA on impedance matching problem 


Lossless impedance matching 
Lossy impedance matching 
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Standard Firefly Algorithm (SFA) 


m Light intensity /(r) is light intensity at a distance r from 
the light source I(r) = he" 
where /, is original light intensity, and y is the light 
absorption coefficient 
m The attractiveness of firefly j to firefly i is defined as 
By = Bo e Ti 
where Bo is the attractiveness at fij = 0 


m The distance between firefly j and firefly / is defined as 
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Standard Firefly Algorithm (SFA) 


2 
xit = Xi + poe A X H G&; 


m First term is a position of firefly í at time t 

m Second term is the exploitation or local search, that is a movement of firefly i 
toward firefly j base on their attractiveness and distance 

m Third term is the exploration or randomness. £; can be specified as a random 
number generated with a value between 0 and 1 





m y and a are very important parameters. 


Ify approaches zero (y > 0), the attractiveness is constant, B,, = Bo, which 
means that every firefly can be seen in any position. 


If the y is nearing infinity or is too large (y + 00), the attractiveness B;, = 0 
and the firefly movements become random 


a is a randomness parameter 
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Modified Firefly Algorithm (MFA) 


m The MFA improves the SFA in three points: 
1) Scaling a computation domain into a unit hypercube 
2) Having a realistic attractiveness between the fireflies 


3) Including a new form of absorption coefficient and a 


randomness parameter 
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MFA — Scaling Domain 


m The computation domain is scaled into a unit hypercube u” 


where 0 < u < 1 ‚and d is the dimension of the searching 
domain. 





Xoia — (Upper Bound — Lower Bound) : X„e„ + Lower Bound 


8 (Xold Lower Bound) 
NEW (Upper Bound—Lower Bound) 
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MFA — Modified Attractiveness 


m The attractiveness of firefly i toward firefly j, Pij is 


modified as > 
Pij = Bolje” 


m Thus, r=1;/Vd 


m where l; a (f(x) = | alla — ma) 
E fmax is the maximum of f(Xj) J=1,..N 
m fmin isthe minimum of f(X) j = 1,.. N 
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MFA — Modified Attractiveness 


m The firefly j is attracted towards the more attractive 
firefly j, the movement is defined as 


x = X; + By - Xt) + Q Ei 
m Thus, 


Xi*" = x? + Polje Y" (X! — XE) + ae 


m The typical value of bg is 1. 
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MFA — New Absorption Coefficient 


m We want the exploration to decrease, but exploitation to 
increase as the number of iterations increase. Thus, the 
parameters y and a should be gradually reduced. Both 
parameters can be formulated as 


Y — Yo: (81)! & = dg: (02)* 


where 0 < 6;, 6> < 1. yg, and agarethe initial value ofy, 
and a „respectively. 
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MFA — New Absorption Coefficient 


m |n practice, we want to have y and a to be very small values at 
the end of iteration. For example, let 


Y — Yo (61) — 1 and Aa = dg: (8,)* = £> 
at t = Max_Iteration, thus, 
= e (log(E1/Y0)/Max_Iteration) „„J 9, = o (log(€2/a9)/Max_Iteration) 


m For example, 
If s; = 0.001, Max_iteration = 50, Yo = 25, we will have 6, = 0.8167 
If €> = 0.001, ag = 0.03 , we will have 6, = 0.8922 
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MFA - Parameters Selection 


m For simplicity, one can use bg = 1; therefore, the two 
parameters to be tuned are yọ > 0 and0 <ap<1. 


m |tcan be seen that y controls the scaling and it indicates how 
fast amplitude can fall, while æ controls the randomness. 

m Consider e-Yor* we define a characteristic length, F, as the 
value of r such that e77/07 = e”1, This means that 


Yo = 1/T*,where0 < T <1 
m Then, if we choose /, so that Yo satisfied the equation above 
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Test Functions Optimization 


m Define the formulas that will be used to compute the 
number of iterations (Max_it) and number of the 
firefly population (N) as follows: 


Max it =20:d +10 
N=5.d+40 


where d is the dimension of problem. 
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Ackley Function 


m The formula of Ackley Function is 


f(x,y) = —20ex p( — 0.24 0.5(x2 + y?) — exp(0.5(cos(27x) + cos(2ny))) +e + 20 


The search space is restricted to —5 < x,y < 5. The global 


minimum of this function is equal to zero, attained at (0,0). 
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Ackley Function 


m The formula of Ackley Function is 


Ackley Function 

















eq 
g 
o» gd 
Y 1 D> Ci s 
a ha s 5 
aei O š 
DA ELO a 
Op Gg 2D {oO 
o <o 
“O gd 3B (© 
DA ARCO 
85 0 5 10 15 20 25 30 35 40 45 50 
iteration 
Locations of fireflies after 10 Final locations of fireflies Comparison of the performance 
iterations after 50 iterations of SFA and MFA for 


minimization of the Ackley 
function with yọ = 25 , and 
a g = 0.02 
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Sphere Function 


m The formula of sphere function is 


f&,y)—x^ ry 


The search space is restricted to —5 < x,y < 5. The global 


minimum of this function is equal to zero, attained at (0,0). 





=f (%, y) 
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Sphere Function 


m The formula of sphere function is 


Sphere Function 





5 5 
4 AŻ 
3h 3y 
2 2 
TH j 1Ę | = 
5 
o 0 z 
m 
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Z -2 
af 8 
-3 R | ` sai EU 
ANN 4 na — 
-5 i E ” ——t e 10° L 4 4 L 4 1 1 1 1 
5 -5 4 3 -2 1 0 1 2 3 4 5 a 5 10 15 20 25 30 36 40 45 50 











iteration 


Locations of fireflies after 10 Final locations of fireflies Comparison of the performance 
iterations after 50 iterations of SFA and MFA with 79 = 150, 
and a@ =0.3 


PENNSTATE 





Rastrigin Function 


m The formula of Rastrigin function is 


f(x,y) = 20 + (x? — 10 - cos(2 + 72x)) + (y* — 10 cos(2 : zy)) 


The search space is restricted to —5 < x,y < 5. The global 
minimum of this function is equal to zero, attained at (0,0). 
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=F (%, y) 
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Rastrigin Function 


m The formula of Rastrigin function is 


4 Rastrigin Function 


“G.G (Gyd 
O 


Best Cost 





L 1 1 L 1 =l L L 4 
0 50 100 150 200 250 300 350 400 450 500 
iteration 


Locations of fireflies after 10 Final locations of fireflies Comparison of the performance 
iterations after 50 iterations, of SFA and MFA with yp = 30, 
converging into (0,0) Ady = 0.1, and Max_it = 500 
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Michalewicz Function 
m Theformula of Michalewicz function is 


| o x? | W^ 
f(x,y) = sin(x); (sin) +sin(y): (sin (y7 


The search space is restricted to —5 < x,y < 5.The global 


minimum of this function is equal to -1.80, attained at 
(2.20,1.57). 
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Michalewicz Function 


m [he formula of Michalewicz function is 


Michalewicz Function 








Best Cost 


























1 L 1 L 1 L L L 1 
a 5 10 15 20 25 30 35 40 45 50 
iteration 


Locations of fireflies after 10 


Final locations of fireflies 
iterations 


Comparison of the performance 
after 50 iterations 


of SFA and MFA with yp = 20, 
and ag = 0.02 
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Outline 


m |ntroduction to impedance matching 
m Nature-inspired optimization algorithms 
m Modified firefly algorithm (MFA) 


m MFA on impedance matching problem 
— Lossless impedance matching 
C Lossy impedance matching 
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Lossless Impedance Matching 


m The Pi3 network is shown below 


Antenna 
Z 





Lossless Pi3 network 


m This network was used for an impedance matching network. The 
load impedance can be obtained from measurement of a 80-meter 
inverted-V dipole antenna between frequencies 3.5-3.85 MHz. In 
simulation the antenna load data is from a record (A15.txt) and a 
source resistance is 50 ohms. 
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Lossless Impedance Matching 
m Transducer Power Gain TPG—-1—|Tr |? 


m J/' isthe reflection coefficient, such that 
l = (Zin — Zs )/( Zin + Zs) 


where Z;, is the input impedance seen by the generator at the equalizer 
input terminal. 


VSWR = a 
Ależ. 


m With y,= 100, a,= 0.6, population = 55, and Max_iteration = 70, an 
average execution time is 30.72 sec. the search space is as follows: 


InF <C, <8nF ,1inF<C,<8nF „and O.l1uH <L < 3 uH 
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Lossless Impedance Matching 
m The optimized component values are obtained as follows: 


C = 4.09 nF, C> = 7.83 nF,and L = 0.6818 uH 


m VSWR of lossless Pi3 network with Min_VSWR = 1.1107, and 
Max_VSWR = 1.3386 


VSWR of Pl network 


—— VSWR 


3.5 


X (MHz.) 
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Outline 


m |ntroduction to impedance matching 
m Nature-inspired optimization algorithms 
m Modified firefly algorithm (MFA) 


m MFA on impedance matching problem 
Lossless impedance matching 
_| Lossy impedance matching 
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Lossy Impedance Matching 





m By adding loss, more bandwidth may be 
acquired at a cost in the transducer power 
gain(TPG) of matching networks. 
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Lossy Impedance Matching - Topology 


lossless lossless 
2-port 2-port 
a) b) 


m Lossy network topology with one resistor 


lossless lossless 
2-port 2-port 
a) b) 


m Lossy network topology with two resistors 


lossless lossless 
2-port 2-port 


m Lossy network topology with three resistors 
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Lossy Impedance Matching 


Antenna 





m [PG can be calculated from 


4Re{Zs}Re{Z, }|Z211° 
|(Zs + 244) (Z, + 222) — 2122211? 


TPG = 


m where Zz is the source impedance, Z, is the load impedance. 214, 212, Z21, 
and Z>» are values of Z-parameter of the matching network. 
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Lossy Impedance Matching 


m Z-parameter of the matching network, that can be obtained from 


1 


Ya = JG6G „A s ESZŁŻ.MA ŻZ= 
cı — JWL1 „Acz To ©) on 


AR = F ak AC, = HM ij AC = F f AL, = a i 


2 "| = AR * AC, * AC * AL, 
he ag y 
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Lossy Impedance Matching 
m We can find a reflection coefficient (T) as following steps 


Za = (jæL : Z,)/(JoL + Z;) 


1 
WN ca 


a 
= Gz? za 5 


ose 
Zin —Z 
P= in S 
Lint Ls 


m The source resistance, Zs = 50 0, and load data are obtained from a 
record (A15.txt). 
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Lossy Impedance Matching 


m T[woobjectives are J; = |I|,andj; = 1 —TPG 
m The final objective function (F) to be minimized, consisting of two 
objectives, is 


F= (Aj, +1 -A)J2) 
Thus, 


F=(AIT|+ (1—2)(1-TPG)) 
where 0 < À < I, in frequency range of 3.5-3.85 MHz. 


m The MFA program run with 79 = 100, ag = 0.6, population = 200, and 
Max_iteration = 60. Each program run has average execution time of 4.30 
minutes. The search space used in the optimization is as follows: 


40 < C, < 2000 pF, 40 < C3 < 2000 pF, 40 < La < 3000 nH,and 1 < R < 1000 Q 
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Lossy Impedance Matching 


m The optimized component values are shown below 


E E 





924.4 1670.7 2681.6 

== 1689.0 1786.2 1930.8 13.8 
M 1397.4 1908.2 2054.0 23.6 
Bzz ei 1199.7 1777.2 2258.2 8.5 
IE 1294.5 1572.3 2237.4 28.1 
ao 886.5 1442.8 2603.1 18.8 
Sas 551.3 1100.9 2824.6 28.0 
1443.0 1509.4 2009.0 28.7 
FZ] 928.8 1285.5 1692.6 40.8 
Sw 1628.5 1503.0 1395.4 42.4 
— 1522.0 1581.5 1098.5 45.9 
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Lossy Impedance Matching 


TPG comparison with different weights 





3.4 3:5 3.6 3.7 3.8 39 4 
Frequency(MHz.) 


TPG of lossy P13 network for A=0, 0.1, ..., 1 
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Lossy Impedance Matching 


WSWR comparison with different weights 





1 M Eee 
3.4 3.5 3.5 af 3.8 3.9 4 
Freguency(MHz.) 


VSWR of lossy P13 network for A=0, 0.1, ..., 1 
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Lossy Impedance Matching 


= Minimum VSWR, maximum VSWR, minimum TPG, and maximum TPG 
result of lossy Pi3 network for weight A = 0,0.1,...,1 


m Min_VSWR Max_VSWR Min_TPG Max_TPG 





1.7101 2.6209 0.9543 0.9600 
1.2693 1.7201 0.9350 0.9446 
Em alas 1.5337 0.9509 0.9894 
Pen 1.1072 1.4748 0.9389 0.9980 
Ea 1.1382 1.4733 0.9112 0.9916 
Em 1.1056 1.4696 0.8494 0.9749 
Pa] 1.1169 1.3468 0.9362 0.9990 
1.0585 1.3006 0.7301 0.9256 
a 1.1006 1.2607 0.7167 0.9197 
| os | 1.0318 1.0908 0.2916 0.4647 
EM 1.0023 1.0080 0.0234 0.0429 
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